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In  previous  papero,  we  have  presented  some  applications  of 
the  principle  of  Invariant  Imbedding  to  radiative  transfer  and 
neutron  diffusion  processes.  This  use  of  Invariance  principles 
was  stimulated  i  y  the  fundamental  v;ork  of  Ambarzumlan  and 
Chandrasekhar,  and  strongly  Influenced  ty  the  point  of  regene¬ 
ration  technique  of  Bellman  and  Harris,  and  the  theory  of 
dynamic  programmj.ng . 

Fundamental  for  the  success  of  these  techniques  as  applied 
to  the  above  processes  Is  the  a'  lllty  to  consider  the  overall 
physical  process  as  a  sequence  of  local  processes.  For  the 
case  of  particles,  this  Is  easily  done.  In  this  paper,  we  wish 
to  Indicate  how  wave  propagation  may  le  considered  In  these 
terns.  It  Is  rather  remarkable  that  our  results  will  le  based 
upon  an  algoi’lthm  ttiat.  In  general,  can  jleld  divergent  series. 

Following  a  provocative  paper  by  Brenmer,  our  aim  Is  to 
show  that  wave  proi.agatlon  can  '  e  discussed  In  terr.s  of  reflec¬ 
tion  and  refraction  at  Infinitesimally  separated  Interfaces. 

We  shall  prove  that  tne  convergence  of  the  Prerrir.er  series  can 
be  established  under  a  simple  assujnptlon  concer’nlng  the  slowly 
varying  nature  of  the  local  wave  num.her. 
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INVARIANT  IMP'i:DDINO,  WAVS  PROPAGATION 
AND  THE  WK3  APPROXIMATION 

Richard  Bellman 
Robert  Kalaba 

1.  Introduction 

In  previous  papers,  we  have  presented  some 

applications  of  the  principle  of  Invariant  1ml  eddlnf^  to 
radiative  transfer  and  neutron  diffusion  processes.  This  use 
of  Invariance  principles  was  stimulated  by  the  fundamental 
work  of  Ambarzumlan  ar^d  Chiandrasekhar,  ,  and  strongly  influ¬ 
enced  by  the  point  of  regeneration  technique  of  Bellman  and 
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Harris,  ,  and  the  theory  of  dynamic  programming, 

Fundamental  for  the  success  of  these  techniques  as  applied 
to  the  above  processes  Is  the  al  lllty  to  consider  the  overall 
physical  process  as  a  sequence  of  local  processes.  For  the 
case  of  particles,  this  Is  easily  done.  In  this  paper,  we  wish 
to  Indicate  how  wave  propagation  may  be  considered  In  tlese 
terms.  It  Is  rather  remarkable  that  our  results  will  be  based 

upon  an  algorithm  that,  In  general,  can  yield  divergent  series. 

0 

Following  a  provocative  paper  by  Bremmer,  '  ,  our  aim  is  to 
show  that  wave  propagation  can  be  discussed  In  terns  of  reflec¬ 
tion  and  refraction  at  Infinitesimally  separated  Interfaces. 

We  shall  prove  that  the  convergence  of  the  Bremmer  series  can 
be  establlo'ied  under  a  simple  assumption  concerning  the  slowly 
varying  nature  of  the  local  wave  num'ter. 

In  addition  to  the  Inslghit  to  the  physical  phenomena 
furnished  by  these  techniques,  we  arc  also  led  to  a  new  metiioJ 
for  studying  the  asymptotic  behavior  of  ordinary  and  partial 
differential  equations  of  linear  type.  For  the  equation 

u"  +  k^( x)u  •  0,  ( 1 • 1 ) 

the  Bremmer  series  Is  an  extension  of  the  known  connection  1 e- 
tween  the  principal  transmitted,  or  refracted,  wave  and  the  WKB 
method.  We  shall  obtain  a  generalization  of  the  WKB  method  to 
vector-matrix  systems  of  the  form 
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y"  +  J{x)y  -  0,  (1.2) 

where  K(x)  Is  a  positive  definite  matrix  for  x  >  0.  This 
yields  a  new  approach  to  the  study  of  the  asymptotic  behavior 
of  the  solutions  of  these  equations, 

Finally,  we  shall  apply  functional  equation  techniques  to 
the  problem  of  determining  the  wave  reflected  l  y  aji  Infinite 
half— plane  of  Inhomogeneous  material.  The  neth.od  we  pursue  Is 
quite  different  from  that  given  In  Lunelerg, 

In  sul  sequent  papers,  we  shall  discuss  the  application  of 
these  Ideas  to  general  hyperbolic  systems,  to  partial  diffe¬ 
rential  equations  of  paratollc  type,  and  to  general  operator 
equations  of  the  foim 

“tt  *  h.j) 

where  A  Is  a  positive  definite  operator.  Furthermore,  we 
shall  consider  spherical  and  cylindrical  geometries. 

2.  A  Localization  Device  for  Wave  Processes 

Consider  a  plane  wave,  ^  ( •  0^  —  wt )  ^  arriving  from  th.e 
homogeneous  space  x  <  0  and  approaching  the  Inhomogeneous 
space  X  >  0.  At  the  boundary,  x  »  0,  the  wave  Is  split 
Into  a  reflected  wave  and  a  refracted  wave.  Let  us  now  suppose 
that  there  Is  an  Immediate  reflected  wave  and  refracted  wave 
obtained  ly  supposing  tl  at  the  Inhomogeneous  medium  Is  actually 
homogeneous  with  wave  nuinl'cr  k(+0).  If  th.e  Inhomogeneous 
region  Is  now  taken  to  ^e  th.e  limit  of  a  sequence  of  Interfaces 
with  this  type  of  reflection  and  refraction  occurring  at  each. 
Interface,  we  can  ol  tain  the  total  reflected  wave  and  the 
disturbance  Inside  the  region  x  >  0  by  adding  up  the  effects 
of  the  reflected  and  refracted  waves  oltalned  In  this  way. 

It  Is  this  principle  which  we  wish  to  verify  under  appro¬ 
priate  assumptions  for  equations  (l.l)  and  (1.2). 

3.  The  Premmer  Series 


Let  us  define  the  sequence 


where 
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MqM  - 


s  l,,^’'k(8)da 


(utr) 


Uj(x)  -  - 


l^/’\(t)clt 

1  /?ook'(s)  /  \  X  .. 

2(k(x))^/^  ''x  k (3)^  ° 


U2^(x) 


1  /7X  k'fs)  f  \  ^3 

Tv^  ?N-a^ 


1  ,/”‘k(t)Jt 
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T7?  i/ 


2(k(x))‘'‘^  “0  k(s) 


ds  , 


l^'’“l:(t)cit 

e  d3. 

g 

11110  oerles  was  derived  by  Bremmer,  ,  using  t!  e  localization 
principle  of  62.  Standard  tecknlqueo  yield 

^  P 

Theorem  1.  ^  lk(x)|  >  a  >0,  for  x  >  0,  and 

/  ®  |k'(x)|dx  <  OD  and  Is  sufficiently  small,  then  tne  Premner 

series ,  u(x)  ■  2  converges  and  represents  a  solution 

n-0  * 

of  (1.1).  Another  linearly  Independent  solution  may  1 e  obtained 
by  replacing  1  ^  -  1 . 


2(k(x)) 


/’OD 


t/ 

X 


k(s) 


(s) 


^2N^ 


4.  Tlie  WKB  Approximation  for  Matrix  Systems 

In  order  to  obtain  an  analogue  of  this  result  for  matrix 
cystems  of  tne  form  of  (1.2),  v;e  introduce  a  matrix  wave— fun 'j t Ion 
(we  s'rall  omit  tlie  scalar  term  )  asroclated  v/ltJ  tl  e 

homogeneous  space  x  <  0  poosesslnr;  t\  c  wave— matrix  K^.  Here 

IKox  10  ^ 

e'*’ denotes  the  matrix  exponential  function;  cf.  .  At  t!  e 

Interface,  x  -  0,  between  x  <  0  and  the  Inhomogeneous  space 
X  >  0,  specified  b\  e  ^  ,  tt  ere  Is  a  reflection  and  a  re¬ 
fraction.  The  dlstur'i  ance  In  x  <  0  Is  given  e^^^^  +  , 

the  original  wave  nlus  a  reflected  wave,  and  there  Is  a  refracted 
wave  In  x  >  0  given  by  where  the  matrices  A  and  P 

are  determined  ly  t!ie  continuity  of  t[ie  functions  and  their 
derivatives  at  x  ■  0.  Hence 
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A 


(Ki  +  Kq)  - 


Kj) 


-  2(K. 


Ko)~^Ko. 


(«.l) 


Taking  account  only  of  the  I'efrac tlonr ,  we  obtain 
Theorem  2 .  Tl'.e  WK^  approximation  to  (1.2)  th.e  solution 
of 

^  -  (lK(x)  -  ^^(x)K'(x))V(x),  V(0)  -  I.  (4.2) 

Under  tie  assumption  tliat  K(x)  Is  positive  definite  for 

X  >  0,  and  that  I  K  '  (x )  |  |  dx  Is  sufficiently  small,  It  can 

0 

be  shown  that  V(x)  Is  the  first  terra  of  a  generalized  Premraer 
series  whlcl  converges  to  a  solution  of  (1.2). 


5.  Reflection  from  an  Inhomogeneous  Space 

Let  us  consider  once  again  the  scalar  plane  wave  case.  We 
wish  to  determine  the  reflected  wave  In  x  <  0  due  to  the 
Inhomogeneous  space  x  >  C.  To  obtain  this,  we  consider  th.e  more 
general  problem  of  determining  the  reflected  wave  from  x  >  z 
due  to  an  Incident  wave  e  O'  '  v'here  ■  k(z  —  0).  Let 

the  coefficient  of  2)  x  <  z  he  u(z).  Tien 

applying  the  localization  principle  of  ^2,  we  see  that 


+  !< ( Z f+'!' ( 1  +  4 ) ' kT , 


which  leads  to  tl  c  Rlccatl  iirferential  equation 

IT  ^  IT 


U'(z)  =  _  1;. 


(-.2) 


If  the  medium  is  ‘'orogeneou.''  for  Zq  <  x ,  with, 
then  u  must  satisfy  tl^  e  ioundary  condition 


u(Zo) 


'■(Zq-O)  - 


0, 
.  i 

0 


k(x) 


("  .3) 


otherwise  we  employ  the^ condition  11m  u(z)  -  0. 

z— ^ 
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